new
restore model.sav    ;ʼģнṹԪڵlink
set echo on
set plot jpg
set plot jpg si 1600 1000
plot set ba wh
plot set cap right
plot set cap si 25
plot set rot 30 0 30

sel group shot1 range group shot_concrete z=(-3.25,0.01)
sel group shot2 range group shot_concrete z=(-4.88,-3.25)
sel group shot3 range group shot_concrete z=(-6.4,-4.88)
sel group shot4 range group shot_concrete z=(-8.94,-6.4)
sel group shot5 range group shot_concrete z=(-10.56,-8.94)
sel group shot6 range group shot_concrete z=(-12.9,-10.56)
sel group shot7 range group shot_concrete z=(-15.00,-12.9)
sel group shot8 range group shot_concrete z=(-17.0,-15.0)
sel group shot9 range group shot_concrete z=(-19.0,-17.0)
sel group shot10 range group shot_concrete z=(-21.0,-19.0)

group excavation range group 12 any group 13 any group 14 any group 15 any
group soil_out range group 16
group exca1 range group excavation z=(-3.25,0.01)
group exca2 range group excavation z=(-4.88,-3.25)
group exca3 range group excavation z=(-6.4,-4.88)
group exca4 range group excavation z=(-8.94,-6.4)
group exca5 range group excavation z=(-10.56,-8.94)
group exca6 range group excavation z=(-12.9,-10.56)
group exca7 range group excavation z=(-15.00,-12.9)
group exca8 range group excavation z=(-17.0,-15.0)
group exca9 range group excavation z=(-19.0,-17.0
group exca10 range group excavation z=(-21.0,-19.0)
;--------------------------------
range name s1 z=(-6.5,0.0)
range name s2 z=(-13.0,-6.5)
range name s3 z=(-17.0,-13.0)
range name s4 z=(-21.0,-17.0)
range name s5 z=(-60.0,-21.0)

;========initial state======================
;model elas
;prop y=40e6 p=0.30  density=1.8e3
model mohr
prop y=40e6 p=0.35 c=8e3 f=10.0 t=1e3 density=1.78e3 range s1
prop y=60e6 p=0.35 c=12e3 f=14.0 t=1e3 density=1.83e3 range s2
prop y=60e6 p=0.35 c=13e3 f=11.0 t=1e3 density=1.76e3 range s3
prop y=60e6 p=0.35 c=21e3 f=21.0 t=1e3 density=1.91e3 range s4
prop y=120e6 p=0.30 c=5e3 f=32.0 t=1e3 density=1.99e3 range s5
;=========================================
fix x range x=(-150.1,-149.9) any x=(149.9,150.1) any
fix y range y=(-0.1 ,0.1 ) any y=(149.9,150.1) any
fix z range z=(-60.1,-59.9)
set gravity=(0.0,0.0,-10.0)
ini szz 0.0 gradient=(0.0,0.0,18e3)
ini sxx 0.0 gradient=(0.0,0.0,9e3)
ini syy 0.0 gradient=(0.0,0.0,9e3)
;=====ʩӵ15kPa====
sel shell apply pres -15e3 range group shell_load
sel shell apply pres -15e3 range group shot1 z=(-0.01,0.01)
;======================================
cyc 1
sel node fix lsys
;-----Ԫڵȷ˳liner, geogrid, pile, cable, shell and beam
;-------------Ƚlink-------------------
sel node ini xpos add 0.0
;========ṹ߽===========
sel node fix y xr zr range y=(-0.01,0.01)
;--------------ṹԪ--------------------------------
sel shell prop iso=(1e-10,0.2) density =0.0 
sel cable prop emod=1e-10 density =0.0 gr_coh=1e20 gr_fri=14.0 gr_k=24e9 &
gr_per=0.314 slide=on slide_tol=1e-6 ycomp=1e20 ytens=1e20
;------------------------------------------
set mech damp combined
sel set damp combined
hist unbal
pl hist 1
step 3000
save initial.sav
;==========ڼ=====================
new
res initial.sav
ini xdisp=0.0 ydisp=0.0 zdisp=0.0
ini xvel=0.0 yvel=0.0 zvel=0.0
sel node ini xdisp=0.0 ydisp=0.0 zdisp=0.0
sel node ini xvel=0.0 yvel=0.0 zvel=0.0
ini state=0
;------------------¶4λ------------------
hist id 10 sel node ydisp id 26950
hist id 11 sel node ydisp id 26934
;--------------A-AһԪ----------------
hist id 101 sel cablesel force cid 178710
hist id 102 sel cablesel force cid 179324
hist id 103 sel cablesel force cid 179872
hist id 104 sel cablesel force cid 180493
hist id 105 sel cablesel force cid 180931
hist id 106 sel cablesel force cid 181312
hist id 107 sel cablesel force cid 181690
hist id 108 sel cablesel force cid 182025
hist id 109 sel cablesel force cid 182331
;========stage1=========== 
model null range group exca1 
step 2000
sel cable prop emod=210e9 density =7.8e3 gr_coh=3.14e4 gr_fri=10.0 gr_k=24e9 gr_per=0.314 &
                         slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6 range group nail1 
sel shell prop iso=(10e9,0.2) density =2.5e3 thick=0.1 ran group shot1 
save stage1.sav
;========stage2=========== 
model null range group exca2
step 2000
sel cable prop emod=210e9 density =7.8e3 gr_coh=3.14e4 gr_fri=10.0 gr_k=24e9 gr_per=0.314 &
                         slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6 range group nail2
sel shell prop iso=(10e9,0.2) density =2.5e3 thick=0.1 ran group shot2
save stage2.sav
;========stage3=========== 
model null range group exca3
step 2000
sel shell prop iso=(10e9,0.2) density =2.5e3 thick=0.1 ran group shot3
save stage3.sav
;========step_excavation=========== 
def step_excavation
loop ii(4,10)
  jj=ii-1
  exca='exca'+string(ii)
  stage='stage'+string(ii)+'.sav'
  nail='nail'+string(jj)
  shot='shot'+string(ii)
  command
    model null range group exca
    step 2000
    sel cable prop emod=210e9 density =7.8e3 gr_coh=3.14e4 gr_fri=10.0 gr_k=24e9 gr_per=0.314 &
                             slide=on slide_tol=1e-6 ycomp=16.5e6 ytens=16.5e6 range group nail 
    sel shell prop iso=(10e9,0.2) density =2.5e3 thick=0.1 ran group shot
    save stage
  endcommand
endloop
end
step_excavation

step 2000
save stage10.sav
